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Abstract 

The notion of density of a finite set is discussed. We proof a general theorem of 
set theory which refines Bose-Einstein distribution. 

In our paper [I] a mistake through negligence occurrecfl- Here we produce a detailed 
proof of the theorem under consideration. 

As in [2], the values of the random variable 3Jl j . . . ? X s £1X6 ordered in absolute value. 
Some of the numbers xi, . . . , x s may coincide. Then these numbers are combined adding 
the corresponding "probabilities", i.e., the ratio of the number of "hits" at x« to the 
general number of trials. The number of equal Xi : Xi = Xi + \ = ■ ■ ■ = Xi + k will be called 
the multiplicity of the value Xi. In our consideration, both the number of trials N and s 
tend to infinity [3]. 

Let Ni be the number of "appearances" of the value x« : X{ < Xi + i, then 

s N 

where M is the mathematical expectation. 

The cumulative probability Vk is the sum of the first k probabilities in the sequence 
Vk = ■h Yli=i Ni, where k < s. We denote NP k = B k . 
If all the variants for which 

s 

J2^ = N (2) 

i=i 

and 

s 

< E, E = MN < Nx, (3) 

i=l 

where x = = q 9iXi , Q = Y^ii=i are equivalent (equiprobable) , then [H[3] the majority of 
the variants will accumulate near the following dependence of the "cumulative probability" 

B<{A r .} = EliA'., 

i=i t=i 
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1 0n page 7 at reference to the paper [2] instead of formula (95) we should refer to formula (94) and 
consider (96), namely v = v' + O(j^jji). This changes the evaluation (4): the exponent should be 3/4 
instead of 1/2. 
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where f3' and v' are determined by the conditions 

B s = N, (5) 

s 

S e/y*, g -?-l = E ' ( 6 ) 

i=l 

as iV — > oo and s — > oo. 

We introduce the notation: .M is the set of all sets {Ni} satisfying conditions (j2J) 
and (jHl); Af{A4} is the number of elements of the set Ai. 

Theorem 1 Suppose that all the variants of sets {N t } satisfying the conditions (G|) 
and (dP are equiprobable. Then the number of variants M of sets {N{\ satisfying condi- 
tions |1P and (T3j) and the additional relation 

iX>-£^3 T ^ iV,3/4+E) < 7 > 

i=i i 

zs less than ClJ ^^ (where C\ and m are any arbitrary numbers, X^=i — an< ^ e ^ s 
arbitrarily small). 

Proof of Theorem 1. 

Let A be a subset of Ai satisfying the condition 



y Ni _ y « — |< A ; 



e f3xi-u _ I 
i=l+l i=l+l 



I I 



V N- V - I < 



A. 



i=i i=i 

where A, (3, v are some real numbers independent of I. 
We denote s 

i=i+l i=l+l 

WNi-Y-^-, — \ = si. 

i=l i=i 

Obviously, if {iVj} is the set of all sets of integers on the whole, then 

s 

N{M \ A} = J2 (®( E ~ E N * X *) S (EU ^(Si - A)Q(S S ^ - A)) , (8) 

{Ni} i=l 

where iVj = iV. 

Here the sum is taken over all integers iVj, 6(x) is the Heaviside function, and 5fc 1: fc 2 
is the Kronecker symbol. 

We use the integral representations 

p —vN fir 

6 NN > = — / ^'Ve^, (9) 

^ J-7T 

1 r°° i 

Qf y ) = — / dx— — e ^ 1+ix \ (10) 
2m J_ OQ x - i 



Now we perform the standard regularization. We replace the first Heaviside function 
in (JED by the continuous function 



for a > 1, y < 

1 _ ePva-ct) for a> i i y > o, 



e /3y(l~a) for a < 0, y < 

1 for a < 0, ?/ > 0, 
where a G (— oo, 0) U (1, oo) is a parameter, and obtain 



1 /"°° 1 1 
©«(y) = 7T- / e"»< 1+to >( : : )dx. (11) 

If a > 1, then Q(y) < Q a {y). 

Let v < 0. We substitute (jUJ) and ffTUj) into (jSj), interchange the integration and 
summation, then pass to the limit as a — > oo and obtain the estimate 

M{M \A}< 



< 



x 



e ~uN+f3E pit 



z'(2tt) 



J -* {N S } 3=1 

Q(S l -A)e(S s ^ l -A)\, J2^ = N, (12) 



where (3 and z/ are real parameters such that the series converges for them. 

To estimate the expression in the right-hand side, we bring the absolute value sign 
inside the integral sign and then inside the sum sign, integrate over ip, and obtain 

-uN+/3E s 

M{M \ A} < 2?r e M~P NiXi + vN $ x 

71 {N t } i=l 

xQ(S l -A)Q(S s ^ l -A). (13) 

We denote 

Z(f3,N) = J2 e ~ PJ:UNlXi ' ( 14 ) 
where the sum is taken over all Ni such that $^i=i Ni = N, 

l s 
i=l i=l+l 

6fofl = (i-^-fe)^ i = l,...,i- 

It follows from the inequality for the hyperbolic cosine cosh(x) = (e x + e~ x )/2 for 
|xi| > 5; \x 2 \ > 5: 

e 5 

cosh(xi) cosh(x 2 ) = cosh(xi + x 2 ) + cosh (a; i — x 2 ) > — (15) 

2 

that the inequality 

l s 

9(S S _* - A)Q(Si — A) < e~ cA cosh (c^iV; - c<p^j cosh(c Ni ~ c ^s-i), (16) 

i=l i=Z+l 



where 

i 



E Hi "T _ \ ^ Hi 

e t3'xi-v' _ l ' 9 s-l— e /3sci-i/ _ ]_ > 

i=l i=Z+l 



holds for all positive c and A. 
We obtain 



M{M \A}< e" cA exp {j3E - vN) 



x 



i i / i \ 

exp{-/? N i x i + v ^2 N i} cosh ( Yl cNi - ^ ) 

{Ni} i=l t=l \i=l / 

s s s 

exp{-/3 A^ + z/ ^ A^j} cosh( cA^ - cfj = 



x 

i=i+l i=Z+l 'i=l+l 

= e^e~ cA x 

x (Ci{u - c, (3) exp(-c0i) + c, (3) exp(c0z)) x 

x {Ca-i(y ~ c, /3) exp(-c0 s _ ; ) + C s -i{v + c, (3) exp(c0 s _,)) . (17) 

Now we use the relations 

d . 

Q^]nCi\p = pt )V=1/ i = <pi] Q^hi( s _i\p = pi iU=u i = <p s _ l (18) 
and the expansion Q{y ± c, (3) by the Taylor formula. There exists a 7 < 1 such that 

ln(C,(f ±c,/3)) = lnC^,^) ±c(lnO)'>^) + ^ (1hCj)> ± jc, (3). 

We substitute this expansion, use formula ([TBI , and see that 0^ is cancelled. 
Another representation of the Taylor formula implies 

In (C,(i/ + c, /?)) = In (0(/3, 1/)) + ~ In ($(/?, 1/)) + 

+ y di/ly + c/P-i/)-^^^*/)). (19) 

A similar expression holds for C, s ~i. 

From the explicit form of the function (i(/3, u), we obtain 

(20) 

where d is given by the formula 

exp(-/3(xi + i/)) 



d 



(exp(-(3(x 1 + 1/)) - 1) 



The same estimate holds for ( s _i. 

Taking into account the fact that CiCs-i — Cs, we obtain the following estimate for 
(3 = (3' and u = v': 

c 2 

M{M \A}< ( s {/3', v') exp(-cA + -l3 2 Qd) exp(£/?' - v'N). (21) 
Now we express Cs(*A fl') i n terms Z(f3, N). To do this, we prove the following lemma. 
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Lemma 1 Under the above assumptions, the asymptotics of the integral 

p —vN rir 

Z(P, N) = — — / dae- lNa ( s ((3 } v + ia) (22) 

has the form 

ZW '^''' \ yi nCMU| (UO| > (23) 
where C is a constant. 

The proof of the lemma readily follows from the saddle-point method and the inequal- 
ities 

Uv + ia,f3)\<Uv,P), \Uv + ia,/3)\<C 8 (vl3), (24) 

which hold, because e u ~ l3xi < 1 for all a ^ 2im, where n is an integer. It follows from 
these inequalities that a = is a saddle point of integral (I2"2"j) [B]. 

We present another proof [7]. 

We have 

Z(J3,N) = / e- iNa ( s {(3,u + ta)da= / e NS ^da, (25) 

27r J_ n 2n J_ n 

where 

s 

S(a,N) = -ia + In ( s ((3,v + ia) = -ia - ^ q { ln[l - e v+ia ~^]. (26) 

8=1 

Here S depends on N, because s, Xi, and v also depend on iV; the latter is chosen so that 
the point a = be a stationary point of the phase S, i.e., from the condition 

We assume that aiiV < s < a2A^, ai,a2 = const, and, in addition, < Xi < B and 
B = const, i — 1, . . . , s. If these conditions are satisfied in some interval f3 £ [0, (3q\ of the 
values of the inverse temperature, then all the derivatives of the phase are bounded, the 
stationary point is nondegenerate, and the real part of the phase outside a neighborhood 
of zero is strictly less than its value at zero minus some positive number. Therefore, 
calculating the asymptotics of the integral, we can replace the interval of integration 
[— 7r, 7r] by the interval [—£,£]. In this integral, we perform the change of variable 

z = y/S(0,N) - S(a,N). (28) 

This function is holomorphic in the disk |a| < e in the complex a-plane and has a 
holomorphic inverse for a sufficiently small e. As a result, we obtain 

E e NS ^da = e NS ^ [ e~ Nz " f{z) dz, (29) 

— e J 7 

where the path 7 in the complex z-plane is obtained from the interval [— e, e] by the 
change (12"8"]1 and 

'dy/S{0,N)-S{a,N) 
da 



(30) 

a=a(z) 



For a small e the path 7 lies completely inside the double sector ie(z 2 ) > c(re z) 2 for some 
c > 0; hence it can be "shifted" to the real axis so that the integral does not change up 
to terms that are exponentially small in N. Thus, with the above accuracy, we have 



e -uN re 



Z(f3,N) = — / e~ Nz f(z)dz. (31) 
2vr J_ £ 

Since the variable z is now real, we can assume that the function f(z) is finite (changing 
it outside the interval of integration), extend the integral to the entire axis (which again 
gives an exponentially small error), and then calculate the asymptotic expansion of the 
integral expanding the integrand in the Taylor series in z with a remainder. This justifies 
that the saddle-point method can be applied to the above integral in our case. 



Lemma 2 The quantity 



1 ^ e -/3E^^ % (32) 



MM) 

V ' {Ni} 

where ^ N{ = N and XiNi < E — i\T 1 / 2+£ ; tends to zero faster than N~ k for any k, e > 0. 

We consider the point of minimum in (3 of the right-hand side of ffTTl) with v(/3, N) 
satisfying the condition 

V - = N 

It is easy to see that it satisfies condition (JSJ). Now we assume that the assumption of the 
lemma is not satisfied. 

Then for £ Nt = JV, x i N i >E- N l ' 2+£ , we have 

{Ni} 

Obviously, [3 provides a minimum of (fT7j) if the assumptions of Lemma 1 are 

satisfied, which contradicts the assumption that the minimum in (3 of the right-hand side 
of ffTTl) is equal to (3' . 

We note that, following the same scheme, G. V. Koval' obtained a direct proof of this 
lemma. 

We set c = N t +a in formula ( |2T|) after the substitution ( |23l) ; then it is easy to see that 
the ratio 

MM \ A) _ 1 
M{M) ~ N m ' 

where m is an arbitrary integer, holds for A = N 3 ^ +£ . The proof of the theorem is 
complete. 

We prove a cumulative formula in which the densities coincide in shape with the Bose- 
Einstein distribution. The difference consists only in that, instead of the set x n of random 
variables or eigenvalues of the Hamiltonian operator, the formula contains some of their 
averages over the cells. In view of our theorem, the which are averages of the energy 
Xk at the ith cell, are nonlinear averages in the sense of Kolmogorov [5] . 

In conclusion, I express my deep thanks to G. V. Koval' and V. E. Nazaikinskii for 
fruitful discussions. 
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